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摘 要 : 人 类 活动 是 土地 利用 变化 的 主要 驱动 因素 ,土地 利用 /有 覆盖 变化 通过 影响 陆地 生态 系统 的 结构 和 功能 间 
接 影 响 生 态 服务 价值 。 本 文 基于 土地 利用 转移 矩阵 和 动态 度 分 析 了 2000 一 2020 年 祁连山 地 区 土地 利用 变化 特征 ， 


应 用 InVEST 模 型 估算 了 祁连山 地 区 生态 系统 的 碳 储量 ,结合 复 利 现 值 公式 估算 了 各 时 期 碳 储 量 经 济 价值 。 结 果 


两 者 之 间 转 换 最 为 显著 。 因 人 类 活动 加 剧 和 城市 化 进程 
HARIN; (2) 由 于 该 地 适宜 的 地 理 环境 ,植被 保护 力度 的 加 强 与 人 工 管理 的 作 / 


表明 :(1) 2000 一 2020 年 草地 和 未 利用 地 是 祁连山 地 区 的 3 


要 土地 利用 类 型 , 占 土 地 利用 类 型 总 面积 80% 以 上 , 且 


的 推进 ,祁连山 地 区 林地 和 未 利用 地 面积 减少 ,其 他 地 类 


] ,2000 一 2020 年 祁连山 地 区 碳 储 


量 增加 44.26x10'*t, 并 呈 "“ 西 北 低 ,东南 高 ”的 空间 分 布 特征 ,与 土地 利用 类 型 分 布 显著 相关 。 其 中 ,耕地 .草地 和 水 
域 碳 储量 增 加 73.1x104t, 林地 和 未 利用 地 碳 储量 减少 28.8x104+t; (3) 2000 一 2020 年 祁连山 地 区 碳 储量 经 济 价值 由 


2325.9x10' 元 增加 至 3908.8x10' 元 , 增 量 为 1382.9x10' 元 ,增长 率 为 68.1%。 其 中 ,草地 碳 储量 经 济 价值 20 a 间 增加 


了 851.8x10' 元 ,增长 率 为 71.5% ,成 为 该 地 区 碳 储量 经 济 价值 增加 的 主要 驱动 力 。 祁 连 山 地 区 碳 储量 呈 增 加 趋势 ， 
未 来 相关 部 门 在 延续 先前 管理 方案 的 同时 ,也 要 警惕 当前 气候 变化 对 碳 储量 的 影响 。 人 研究 结果 可 为 区 域 制定 优化 


土地 利用 目标 推动 可 持续 发 展 和 应 对 全 球 气候 变化 提供 科学 依据 。 


关键 词 : 土地 利用 ; InVEST 模 型 ， 碳 储量 ; 碳 储量 经 济 价值 ;祁连山 


全 球 气候 变 暧 引发 的 一 系列 环境 问题 已 对 生 
态 环 境 和 社会 经 济 的 可 持续 发 展 产 生 深 远 影响 中。 
人 类 活动 造成 的 土地 利用 变化 / 履 盖 导致 陆地 生态 
系统 碳 源 / 汇 失 衡 成 为 大 气温 室 气 体 浓 度 增 加 的 主 
要 原因 之 一 。 土 地 利用 变化 通过 改变 陆地 生态 系 
统 结构 和 功能 ,影响 陆地 生态 系统 的 碳 储量 ” ,间接 
影响 碳 储量 经 济 价 值 。 因 此 ,研究 区 域 土地 利用 变 
化 对 碳 储存 与 碳 储量 经 济 价值 的 影响 ,可 为 我 国 制 
定 减 排 增 汇 政策 提供 重要 的 科学 基础 ,对 实现 “ 双 
矶 "目标 具有 重要 意义 。 

当前 已 有 大 量 学 者 对 土地 利用 变化 和 碳 储量 
的 估算 进行 了 研究 ,传统 的 碳 储量 佑 算 方法 主要 有 
样 地 清查 法 .野外 获取 法 和 IPCC 清单 法 等 ? joc we 
方法 获取 数据 路 径 可 靠 ,估算 结果 精确 ,但 因 普 适 
PED .采样 周期 长 和 工作 繁琐 等 不 足 ,无 法 进行 
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尺度 研究 中。 随 着 技术 的 发 展 ,多 数学 者 基于 Cen- 
tury 模 型 .RothC 模型 和 LPJ-GUESS 模型 等 5 9 进行 
区 域 碳 储量 估 算 。 虽 然 模 型 估算 方法 弥补 了 传统 
方法 的 不 足 , 但 多 数 集中 于 单一 土地 利用 类 型 的 研 
究 , 忽 略 了 各 土地 利用 类 型 之 间 存 在 相互 竞争 的 复 
杂 关 系 ""。 然 而 ,InVEST 模 型 不 仅 可 以 弥补 传统 人 
算 方 法 的 不 足 , 而 且 具 有 参数 少 、 模 型 精度 高 和 运 
行 速度 快 等 优点 ,已 被 国内 外 大 量 学 者 应 用 于 不 
同 尺度 综合 土地 利用 变化 与 碳 储 量 之 间 复 杂 的 响 
应 机 制 研 究 中 。 如 李 俊 等 中 基于 InVEST 模 型 对 
昆明 市 的 碳 储量 进行 评估 ,研究 表明 20 a 间 因 建筑 用 
地 的 扩张 导致 碳 储量 下 降 了 9x10't, 为 昆明 市 优化 土 
地 利用 和 可 持续 发 展 目标 提供 数据 支撑 ;Li 等 中 基于 
InVEST 模 型 对 长 春 市 土地 利用 与 碳 储 量 的 变化 进 
行 研究 ,结果 表明 , 近 10a 间 ,耕地 向 建设 用 地 转移 
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是 该 地 区 碳 储 量 损失 的 主要 原因 。 通 过 测算 不 同 
土地 利用 方式 下 的 碳 储量 变化 ,并 将 其 转化 为 碳 排 
放 减 少 或 碳 吸 收 增加 对 应 的 经 济 价值 ,可 以 对 土地 
利用 变化 的 效益 进行 量化 评估 ,这 种 估算 的 结果 通 
常 作 为 政策 制定 和 资源 配置 的 决策 参考 。 例 如 , 伍 
格致 等 引 利 用 造林 成 本 法 和 碳 税 法 的 算数 平均 值 
对 湖南 省 森林 的 碳 储量 经 济 价值 进行 估算 ,为 湖南 
省 生态 保护 机 制 提 供 数据 基础 。 值 得 注意 的 是 , 国 
内 碳 交 易 市 场 不 成 熟 ,尚未 形成 统一 的 碳 储量 经 济 
价值 估算 方法 ,主要 有 市 场 价值 法 .造林 成 本 法 、 人 
工 固定 CO; 成 本 法 、 碳 税法 .支付 意愿 法 和 实物 期 权 
法 等 中 ,这 些 方法 理论 通俗 易 懂 、 所 需 数据 可 通过 
调查 .统计 年 鉴 和 相关 资料 获得 ,但 通常 忽略 了 货 
币 的 时 间 价 值 ,往往 通过 某 一 时 期 或 固定 的 碳 交 易 
价格 去 估算 不 同时 期 的 碳 储 量 经 济 价值 ,从 而 导致 
估算 结果 精确 度 较 差 。 碳 税法 采用 瑞典 碳 税 率 ， 
但 瑞典 碳 税率 对 于 我 国 偏 高 ,不 符合 我 国 现 发 展 阶 
段 ,会 使 估算 的 结果 偏 高 “。 因 此 ,在 佑 算 不 同时 
期 的 碳 储量 经 济 价值 时 , 碳 价格 的 时 间 价 值 非常 重 
要 。 目前 针对 这 个 问题 ,国外 已 有 相关 人 研究"。 
例如 ,Hussainzad 等 "利用 净 现 值 法 ,评估 绢 享 州 固 
碳 的 经 济 价值 ,为 可 持续 环境 保护 和 低 碳 经 济 转型 
提供 有 力 的 支持 。 然 而 ,国内 关于 该 问题 的 研究 
相对 欠缺 。 因 此 ,本 研究 采用 财务 管理 .金融 学 和 
工程 经 济 学 领域 广泛 使 用 的 复 利 现 值 法 ,通过 中 
央 银 行规 定 的 贴现 率 ,考虑 了 碳 价格 变 化 的 时 间 
MMES ,对 祁连山 地 区 不 同时 期 碳 储 量 的 经 济 价值 
进行 估算 。 

祁连山 地 区 是 中 国 西北 地 区 重要 的 生态 安全 
屏障 ,在 维系 区 域 生态 平衡 中 扮演 着 不 可 或 缺 的 角 
色 呈 。 但 随 着 近年 来 人 口 的 快速 增长 .经济 的 迅速 
发 展 以 及 城镇 化 进程 的 急速 推进 ,该 地 区 土地 利用 
变化 明显 ,一 定 程度 上 改变 了 人 研究 区 的 碳 储 存 功 
能 。 现 有 的 研究 主要 关注 不 同 影响 因子 与 生态 环 
境 系统 之 间 的 响应 机 制 “” ,通过 量化 生态 系统 服 
务 价值 ,可 以 更 好 的 建立 生态 补偿 标准 ,并 有 效 权 
衡 生态 服 务 供需 平衡 效应 2 ,而 对 生态 系统 碳 储量 
经 济 价值 的 估算 仍 未 知 。 基 于 此 ,本 研究 使 用 
2000 一 2020 年 祁连山 地 区 的 五 期 土地 利用 数据 , 通 
过 土地 利用 转移 矩阵 和 动态 度 分 析 了 土地 利用 变 
化 特征 ,应 用 InVEST 模型 评估 了 该 地 区 2000 一 
2020 年 陆地 生态 系统 碳 储量 ,并 结合 中 央 人 民 银 行 


的 贴现 率 ,使 用 复 利 现 值 公式 估算 了 碳 储量 经 济 价 
值 。 本 研究 揭示 了 土地 利用 变化 如 何 影 响 碳 储量 
并 佑 算 了 碳 储量 经 济 价 值 , 对 于 制定 土地 利用 策 
略 、 优 化 生态 环境 .应 对 气候 变化 及 推动 碳 交 易 市 
场 发 展 具有 重要 现实 意义 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

ASI LL Hh X (93.4° ~103.4° E, 35.8°~40.0°N) 位 
于 中 国 西北 部 ,地 处 甘肃 .青海 省 交界 (图 1) ,总 面 
积 约 为 19.1x10! hn, 东 起 乌 鞘 岭 , 西 止 当 金山 口 , 北 
临河 西 走 廊 , 南 依 某 达 木 岔 地。 该 地 具有 典型 大 陆 
性 气候 特征 ,是 我 国生 物 多 样 性 保护 热点 区 域 和 重 
要 的 生态 环境 区 域 ”。 祁 连 山 地 区 土壤 类 型 多 样 , 具 
有 明显 的 垂直 分 布 特征 , 且 矿 产 和 人 文 资 源 丰 富 汪 。 
因此 ,祁连山 的 生态 保护 不 仅 关 乎 地 区 生态 的 稳 
定 , 更 涉及 我 国 整 个 西部 地 区 的 生态 平衡 与 经 济 社 
会 的 可 持续 发 展 。 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2020)4619 
号 ,对 底 图 无 修改 。 
图 1 祁连山 地 区 示意 图 
Fig. 1 Map of Qilian Mountains 


1.2 数据 来 源 

土地 利用 数据 来 源 于 中 国 科学 院 资 源 环境 科 
学 与 数据 中 心 (https://www.resdc.cn/) ;全 国 年 平均 气 
候 数据 源 于 《2020 中国 气候 公报 》; 人 研究 区 的 年 平均 
气温 数据 来 自 4 中 国 地 面 资 料 气候 标准 值 数据 集 》 
(http://data.cma.cn) 。 碳 密度 数据 如 表 1 所 示 ,主要 
基于 两 个 来 源 : 第 一 ,依据 参考 文献 求 得 实测 数据 
的 平均 值 ”, 第 二 ,针对 难以 获得 的 数据 ,利用 校正 


11 期 雷 ” 世 等 :祁连山 地 区 土地 利用 


表 1 祁连山 地 区 不 同 土地 利用 类 型 碳 密 度 
Tab.1 Carbon density of different land types 


in Qilian Mountains 


土地 利用 碳 密度 /(t: hm) 

类 型 ” 地 上 生物 量 地 下 生物 量 土壤 死亡 有 机 质 总 碳 密度 
耕地 4.56 7.45 106.30 3.67 121.98 
林地 51.08 37.51 174.95 6.71 270.25 
草地 1.02 6.15 128.51 3.30 138.98 
水 域 0.00 0.00 170.27 0.00 170.27 
建设 用 地 0.00 0.00 0.00 0.00 0.00 
未 利用 地 0.12 1.57 52.96 0.00 54.65 
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据 进行 修正 而 获得 。 由 于 建筑 用 地 土地 硬化 ,本 研 
究 将 建筑 用 地 的 碳 密度 视 为 0。 

1.3 研究 方法 

1.3.1 土地 利用 转移 矩阵 ”土地 利用 转移 矩阵 是 在 
一 定时 间 和 区 域内 土地 利用 类 型 相互 转换 的 二 维 矩 
阵 , 表 示 土 地 利用 结构 和 土地 利用 类 型 的 变化 方向 ””。 


Sa Si, 
a ae (1) 
Sat Sia Spn 


式 中 :5; 为 某 时 期 第 i 类 土地 利用 类 型 转化 为 第 j 类 
土地 的 面积 (hm?);n 为 土地 利用 类 型 数量 。 
13.2 土地 利用 动态 度 ”土地 利用 动态 度 是 反应 区 
域 一 定时 间 内 土地 利用 变化 速率 的 指标 ,包括 单一 
土地 利用 动态 度 (K) 和 综合 土地 利用 动态 度 (LC)"”。 
U,- 

U, 
式 中 :kK 为 某 类 土地 利用 类 型 动态 度 ;U6, 和 U0 分别 为 
研究 初期 和 末期 某 土 地 利用 类 型 的 面积 (hm );7 为 
研究 时 长 (年 )。 

YAW, 


LC =|=——_— 

2J LU, 
AP LC 为 土地 利用 综合 动态 度 ;7U 为 研究 起 始 时 
间 第 ;类 土地 利用 类 型 的 面积 (hm ) ;ALU 为 研究 期 
内 第 i 类 土地 转化 为 第 j 类 土地 利用 类 型 面积 的 绝 


y= 


K= 


a x $ x 100% (2) 


xL x100% (3) 
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对 值 ;7 为 研究 时 长 (年 )。 
1.3.3 InVEST 模 型 ”本 研究 采用 InVEST 模 型 中 的 
Carbon storage and sequestration 模块 进行 碳 储 量 计 
算 ,陆地 碳 库 主要 包括 :地 上 生物 量 、 地 下 生物 量 、 
土壤 和 死亡 有 机 质 四 大 碳 库 "”。 碳 储量 由 各 土地 
利用 类 型 碳 密度 与 其 相应 的 土地 利用 类 型 面积 相 
乘 所 得 。 计 算 方法 如 下 所 示 : 

Cp = C atore Crete t+ Co * Caead (4) 


n 
Coui = yA, x Ci 
k=1 


式 中 ;Gi 为 某 土 地 利用 类 型 碳 储 量 之 和 (t) ; C 为 
地 上 生物 量 碳 密度 (thm”); Coo 为 地 下 生物 量 碳 
密度 (thm-2); C, NERE CG- hm’); C N 
死亡 有 机 质 碳 密度 (thm”) ;4; 为 某 土地 利用 类 型 面 
FR Chm?) ;Ci 为 某 土地 利用 类 型 磋 密 度 (t+hm?) ;Ca 
为 区 域 陆地 生态 系统 总 碳 储量 (t)。 
13.4 复 利 现 值 法 复 利 现 值 法 产生 于 商品 经 济 学 
中 ,涵盖 了 货币 的 时 间 价 值 ,是 个 人 和 企业 进行 投 
资 决 策 时 的 重要 指标 。 尽 管 该 方法 在 经 济 学 中 被 
广泛 应 用 ,但 在 陆地 碳 储量 经 济 价值 估算 方面 还 未 
被 探索 。 本 研究 通过 碳 交 易 平台 获得 2019 一 2021 
年 平均 碳 交 易 价 格 作为 2020 年 的 碳 价格 ,通过 复 利 
现 值 公式 和 转换 系数 得 到 2000 年 .2005 年 .2010 年 
和 2015 年 的 碳 交易 价格 ,如 表 2 所 示 。 

复 利 现 值 表 示 未 来 一 定时 间 的 特定 
利 计算 的 现 值 。 

P,=Fx (P/F,e,n) (6) 

式 中 :PP, 为 第 n 期 的 碳 价 格 (元 ); 了 为 碳 价格 终 值 
(元 ); (P/F,eyn) 为 复 利 现 值 系 数 ;e 为 贴现 率 ;n 为 
期 数 。 
13.5 碳 储 量 价 值 估 算 ” 碳 储量 价值 表示 生态 系统 
吸收 并 储存 碳 所 需要 的 经 济 价值 。 


P= $C,” P, (n=1,2,...,n) (7) 


式 中 :P. 为 陆地 生态 系统 总 碳 储 量 价 值 (元 );C, 为 第 
n 期 矶 储量 (t) ;P, 为 第 n 期 碳 交 易 价 格 (元 中)。 


(k=1,2,.…,6) (5) 


资金 按 复 


表 2 2000 一 2020 年 碳 交易 价格 
Tab.2 Carbon price from 2000 to 2020 


年 份 2000 年 2005 年 


20104F 2015 年 2020 Œ 


碳 交 易 价 格 /元 上) 105 116 


128 156 173 
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2 结果 与 分 析 


土地 利用 时 空 变 化 分 析 

通过 祁连山 地 区 土地 利用 面积 得 到 祁连山 地 
区 2000 一 2020 年 土地 利用 转移 矩阵 和 动态 度 变化 
(图 2, 表 3) ,祁连山 地 区 主要 以 草地 和 未 利用 地 为 
主 , 面 积 占 比 约 83% ,各 土地 利用 类 型 间 转 换 频 
繁 。 草 地 面积 增加 最 显著 ,在 过 去 20 a 增加 了 
3347.5 hm ,增长 率 为 4.1% ,其 中 ,第 一 个 10 a 草地 
面积 增加 了 3605.2 hm ,增幅 为 4.1% ,而 第 二 个 10 a 
呈现 减少 趋势 ,减少 了 257.7 hwo IKER, 耕地 和 建 
筑 用 地 ,在 过 去 20 a 分 别 增加 1386.4 hm’? 241.2 hm? 
和 136.1 hm2 ,增幅 分 别 为 17.8% .4.3% 和 32.4%。 然 
而 ,林地 和 未 利用 地 面积 在 2000 一 2020 呈现 减少 趋 
势 ,其 中 ,未 利用 地 面积 减少 明显 ,减少 了 5070.5 
个 10 a 减少 了 4734.7 hm ,降幅 为 6% ,第 


2.1 


A 
hm ,第 一 


(a) 2000 一 2010 年 /hm? 


81600 70000 


(b) 2010 一 2020 年 /hm? 


00 5730 


T> 


85200 30000 70000 


二 个 10 a 减少 了 335.8 hm ,降幅 为 0.5%。 林 地 在 过 
ZF 20 a 共 减少 40.7 hm ,降幅 最 低 ,为 0.2%。 根 据 
2000—2020 年 单一 土地 利用 动态 度 ( 图 3) ,建筑 用 
地 动态 度 绝对 值 最 大 ,为 3.24%, 其 次 为 水 域 .未 利 
用 地 和 耕地 。 然 而 ,林地 动态 度 绝对 值 最 小 , 仅 为 
0.02% 。 从 综合 土地 动态 度 来 看 ,祁连山 地 区 
2000 一 2020 年 土地 动态 度 呈 波动 式 变 化 ,由 2000 年 
的 0% 增 加 到 2010 年 的 0.25% ,随后 下 降 到 2015 年 
的 0.01% ,最 后 在 2020 年 回升 至 0.03% ,整体 呈 增 加 
趋势 。 

祁连山 地 区 土地 利用 间 分 布 呈 现 不 同 
的 格局 (图 3) ,耕地 面积 占 比 约 2.98% , 主要 分 布 在 
东北 部 和 东南 部 ,变化 区 域 主要 分 布 在 东南 部 ; 草 
地 面积 占 比 40% ,主要 分 布 在 北部 .南部 和 东南 部 ， 
变化 区 域 主要 分 布 在 北部 和 东南 部 中 心地 带 ; 林 地 
在 过 去 20a 有 0.2% 的 小 幅度 减少 ,主要 分 布 在 东南 
部 ,变化 区 域 零 散 分 布 在 研究 区 ;水 域 主要 分 布 在 


类 型 空 


(c) 2000 一 2020 年 /hm2 


0 79200 5560 


7000! 


OS 耕地 国 国 林地 Deh Bh 国 国 建筑 用 地 未 利用 地 
图 2 祁连山 地 区 土地 利用 转移 矩阵 
Fig. 2 Land use transfer matrix in Qilian Mountains 
表 3 祁连山 地 区 土地 利用 类 型 变化 
Tab.3 Land use type change in Qilian Mountains 
TER 2000—2005 4: 2005—2010 4 2010—2015 Æ 2015—2020 4 2000—2020 4 
pe 变化 eaen 变化 变化 - 变化  - 变化 . 
类 型 bag 态度 态度 7 态度 = 动态 度 态度 
幅度 /hm 动态 度 /% 。 幅度 /hmz EE 幅度 /hm POEA gme 动态 度 /% g 动态 度 /% 
耕地 60.1 0.11 107.6 0.19 -11.2 -0.02 84.8 0.15 241.2 0.43 
林地 -6.3 0.00 -50.8 -0.03 36.6 0.02 -20.2 -0.01 -40.7 -0.02 
草地 -69.8 -0.01 367.5 0.45 27.5 0 -28.5 -0.03 334.7 0.41 
水 域 6.0 0.01 953.9 1.22 -59.0 -0.07 485.5 0.56 138.6 1.78 
建筑 用 地 2.8 0.07 56.2 1.33 66.5 1.39 10.6 0.19 136.1 3.24 
未 利用 地 7.1 0.00 -47.4 -0.6 -60.5 -0.01 -27.5 -0.04 -5070.5 -0.64 
综合 动 0.00 0.25 0.01 0.03 0.27 
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变化 对 碳 储量 的 影响 及 经 济 价值 估算 


“耕地 2.91% ”硬是 水 域 4.09% 
国王 林地 8.66% 。。 国 国 建筑 用 地 0.22% 
“草地 42.7% 未 利用 地 41.42% 
(c) 2020 年 


“” “耕地 3.04% 
园林 地 8.64% 
草地 44.45% 


E 7K R4.81% 
E 建筑 用 地 0.29% 
未 利用 地 38.77% 


”耕地 3.00% 
ME 林地 8.63% 
m 草地 44.59% 未 利用 地 38.95% 
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(b) 2010 年 


À 


0 


EE 水 域 4.59% 
EE 建筑 用 地 0.25% 


(d) 2000 一 2020 年 土地 利用 变化 


Zz 


0 200 km 


不 变 区 域 E 林地 变化 区 域 E 水 域 变化 区 域 
未 利用 地 变化 区 域 ““” 耕地 变化 区 域 MN 草地 变化 区 域 


E 建筑 用 地 变化 区 域 


图 3 祁连山 地 区 土地 利用 空间 分 布 


Fig. 3 Spatial distribution map of land use in Qilian Mountains 


南部 和 西北 部 中 心 ,变化 区 域 零 散 分 布 于 西北 部 ， 
小 部 分 集中 分 布 于 东北 部 的 中 段 ;截至 2020 年 , 建 
筑 用 地 面积 增加 ,增加 区 域 主要 分 布 于 研究 区 北部 
和 中 段 ,零散 分 布 于 西部 和 东南 部 。 
2.2 碳 储量 时 空 变 化 分 析 

由 祁连山 地 区 各 土地 利用 类 型 碳 储量 的 时 间 
变化 趋势 可 知 (图 4) ,祁连山 地 区 总 碳 储量 呈 波 动 
上 升 趋势 ,从 2000 年 的 22.2x105t 增 加 到 2020 年 的 
22.6x10t, 增 幅 为 2%。 值 得 注意 的 是 ,在 2000 一 
2005 年 和 2010 一 2015 年 总 碳 储量 有 小 幅度 的 下 
降 。20 a 间 耕地 .草地 和 水 域 碳 储量 增加 ,林地 和 未 
利用 地 碳 储量 减少 。2000 一 2020 年 耕地 碳 储量 呈 
增加 一 减少 一 增加 的 变化 趋势 ,20 a 间 增加 了 29.4x 
10:t ,增幅 为 4.3% , FEP , 2010—2015 年 碳 储量 下 降 
了 1.4x10:t ,降幅 为 0.2%。 林 地 碳 储量 呈 减 少 一 增 
加 一 减少 的 变化 趋势 ,20 a 间 碳 储量 下 降 了 11x10; 
t, 降幅 为 0.3% ,其 中 ,2010 一 2015 年 碳 储 量 增加 
9.9x10° t, 增 幅 为 0.2%。 草 地 碳 储 量 呈 减少 一 增加 
一 减少 的 变化 趋势 ,20 a 间 增加 了 465.2x10:t ,增幅 
为 4.1%, 其 中 ,2000 一 2005 年 和 2015 一 2020 年 碳 储 


量 分 别 下 降 97.0x10’t 和 39.6x10; ,降幅 分 别 为 
0.1% 和 0.3%。 水 域 碳 储 量 呈 增加 一 减少 一 增加 的 
变化 趋势 ,过 去 20 a 增加 了 236.1x10’ t, 1% 1H 
17.8% , 其 中 , 2010—2015 年 碳 储量 下 降 ,降幅 为 
0.7% 。 未 利用 地 碳 储量 20 a 间 下 降 了 227.1x10; t, 
降幅 为 6.4%。 

从 祁连山 地 区 2000 一 2020 年 各 土地 利用 类 型 
碳 储量 空间 分 布 情况 可 知 (图 5)。 不 同时 期 祁连山 
地 区 碳 储量 空间 分 布 大 致 相同 ,总 体 呈 现 出 "西北 
低 东 南 高 ”的 空间 分 布 格局 。 碳 储量 高 值 区 主要 分 
布 在 东部 和 中 部 ,以 林地 和 草地 为 主 。 碳 储量 中 值 
区 零散 分 布 在 西部 ,以 草地 和 水 域 为 主 。 碳 储量 低 
值 区 主要 分 布 在 西部 和 中 部 ,以 未 利用 地 为 主 。 
2000—2020 年 祁连山 地 区 总 碳 储量 增加 区 域 主 要 
分 布 在 中 部 区 域 , 碳 储量 减少 区 域 零星 分 布 在 祁 连 
山 整个 地 区 。 
2.3 碳 储量 经 济 价值 变化 

祁连山 地 区 碳 储量 经 济 价值 从 2000 年 2325.9x 
10° 70 34 K Z 2020 4 HY) 3908.8 x 10°76 , HK XH 
68.1%( 表 4)。 耕 地 在 过 去 20 a 碳 储量 经 济 价值 增 
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图 4 祁连山 地 区 2000 一 2020 年 各 土地 利用 类 型 碳 储量 变化 
Fig. 4 Change of carbon storages by land use type from 2000 to 2020 in Qilian Mountains 
(a) 2000 年 N (b) 20104 À 
碳 储量 /t 
270.25 一 
0 0 200 km 0 200 km 
S| | 
(c) 2020 年 人 (d) 2000 一 2020 年 碳 储量 变化 À 
碳 储量 ( = Zk 碳 储量 /t 
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| L 
图 5 祁连山 地 区 碳 储量 空间 分 布 
Fig. 5 Spatial distribution of carbon storages in Qilian Mountains 
表 4 祁连山 地 区 2000 一 2020 年 各 土地 利用 类 型 碳 储量 经 济 价值 
Tab. 4 Value of carbon storages in various regions from 2000 to 2020 in Qilian Mountains /10' 元 
土地 利用 类 型 2000 年 2005 年 2010 年 2015 年 2020 年 
耕地 71.3 79.6 98.5 108.8 122.5 
林地 469.7 518.8 570.7 697.0 772.0 
草地 1191.0 1314.7 1516.0 1848.3 2042.8 
水 域 139.6 154.4 191.1 231.4 270.9 
未 利用 地 454.3 502.0 520.7 634.1 700.6 
总 价值 2325.9 25693.0 2888.0 3519.6 3908.8 
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加 51.2x10 元 ,增幅 为 71.8% , EP , 2010—2015 年 
增幅 最 大 ,为 21.7%。 林 地 和 草地 碳 储量 经 济 价值 
分 别 增加 302.3x10' 元 和 851.8x10 元 ,增幅 分 别 为 
64% 和 71.2% ,同样 在 2010 一 2015 年 增幅 最 高 ,分别 
为 22.2% 和 21.9%。 水 域 碳 储 量 经 济 价值 在 20 a [ia] 
增加 131.3x10* 元 ,增幅 为 94% ,但 在 2000 一 2005 年 
增幅 呈 最 高 ,为 23.8%。 未 利用 地 碳 储量 经 济 价值 
增加 246.3x10 元 ,增幅 为 53.9%。 根 据 各 土地 利用 
类 型 碳 储 量 经 济 价 值 占 比 (图 6), 草 地 、 林 地 和 未 利 
用 地 碳 经 济 价 值 占 比 最 高 ,平均 占 比分 别 为 51.9%、 
19.9% 和 18.6%。 相 反 ,耕地 碳 储量 经 济 价值 占 比 最 
小 ,年 平均 占 比 为 3.1%。 


国 未 利用 地 Oka 国 草 地 国 林 地 国 耕地 


碳 储量 经 济 价值 占 比 /% 
二 


2000 2005 2010 2015 2020 
年 份 


图 6 祁连山 地 区 2000 一 2020 年 各 土地 利用 类 型 碳 储量 
经 济 价值 占 比 
Fig. 6 Percentage distribution of carbon economic from 
2000 to 2020 in Qilian Mountains 


3 讨论 


3.1 碳 储量 对 土地 利用 变化 的 响应 

长 期 以 来 ,土地 利用 变化 是 影响 陆地 生态 系统 
中 碳 储存 变化 的 关键 因素 ,无 论 是 林地 到 耕地 的 转 
变 , 还 是 建筑 用 地 到 农业 用 地 的 转型 , 亦 或 是 未 利 
用 地 的 开发 和 建设 ,这 些 变化 对 陆地 生态 系统 碳 存 
储 和 时 空 变 化 产生 深远 的 影响 中。 本 研究 明确 了 
祁连山 地 区 2000 一 2020 年 土地 利用 变化 特征 ,结合 
InVEST 模 型 的 Carbon storage and sequestration 模块 
分 析 了 祁连山 地 区 各 土地 利用 类 型 碳 储量 时 空 分 
布 特征 。 随 着 人 口 的 增加 和 城镇 化 的 推进 ,祁连山 
地 区 的 土地 利用 模式 发 生 较 大 的 转变 。 了 解 土地 
利用 变化 对 碳 储量 的 影响 ,有 助 于 制定 合理 的 土地 
管理 政策 ,为 减缓 气候 变化 提供 数据 支撑 ”。 
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众所周知 ,耕地 和 林地 土壤 有 机 碳 密度 高 , 具 
有 和 较 强 的 碳 储 存 能 力 ,而 建筑 用 地 土壤 碳 密度 低 ， 
相应 的 碳 储 存 能 力 较 弱 。 土 地 利用 变化 通过 改变 
生态 系统 分 布 格局 和 土壤 呼吸 速率 等 影响 陆地 生 
态 系统 的 碳 储量 。 例 如 , 李 瑾 正 等 ”基于 InVEST 模 
型 估算 湖北 省 陆地 生态 系统 碳 储 量 ,研究 表明 ， 
1990—2015 年 主要 因 林 地 等 碳 密 度 值 高 的 地 类 转 
化 为 建筑 用 地 导致 碳 储 量 减 少 了 29.23 Tg, 说 明 建 
筑 用 地 的 扩张 规模 直接 影响 碳 储量 变化 速率 ;朱文 
博 等 基于 土地 利用 变化 情景 的 生态 系统 碳 储量 
评估 ,研究 表明 ,2005 一 2015 年 建设 用 地 的 肆意 扩 
张 和 林 地 的 退化 是 导致 碳 密 度 下 降 的 主要 原因 ,对 
区 域 土地 利用 的 合理 规划 具有 重要 意义 。 本 研究 
中 草地 作为 祁连山 地 区 主要 的 土地 利用 类 型 并 成 
为 主要 的 碳 库 ,对 陆地 生态 系统 起 到 重要 作用 。 与 
邓 蔡 等 对 祁连山 地 区 碳 储 量变 化 的 研究 得 出 一 
致 的 结论 , 即 祁连山 地 区 碳 储量 整体 呈 增 加 趋势 ， 
这 主要 归 因 于 两 方面 ,第 一 , 碳 密度 值 低 的 未 利用 
地 转化 为 碳 密 度 值 高 的 草地 ,水域 和 耕地 等 地 类 ， 
增加 了 陆地 生态 系统 碳 储 存 能 力 ”。 第 二 ,国家 颁 
布 的 野生 动物 保护 法 和 实施 天 保 工 程 所 采取 的 封 
山 育 林 、 禁 止 放牧 等 措施 ,减少 人 类 活动 ,使 固 碳 能 
力 较 强 的 土地 在 一 定 程度 上 得 到 了 保护 。 然 而 ,在 
2010 一 2015 年 祁连山 地 区 碳 储量 出 现 小 幅度 下 降 ， 
这 是 因为 城镇 化 速度 加 快 ,建筑 用 地 扩张 占据 了 大 
量 的 水 域 和 未 利用 地 ,导致 该 期 间 碳 储量 下 降 。 由 
此 可 见 , 其 他 地 类 转化 为 建筑 用 地 对 陆地 碳 储 存 能 
力 的 降低 程度 影响 较 大 。 

另外 ,中 国 西北 祁连山 地 区 作为 气候 变化 敏感 
区 域 , 其 生态 环境 脆弱 ,受到 气候 变化 和 人 类 活动 
的 双重 影响 。 尽 管 本 研究 表明 当地 生态 系统 碳 储 
HELE 2000—2020 年 整体 呈 上 升 趋势 ,但 是 由 于 特殊 
的 地 理 特征 和 气候 条 件 , 该 地 区 往往 面临 土地 退化 
和 生物 多 样 性 丧失 等 问题 。 为 此 ,建立 长 期 的 土 
地 利用 变化 监测 体系 十 分 必要 。 随 着 未 来 气候 变 
化 的 加 剧 和 经 济 社会 发 展 的 需求 ,土地 利用 变化 程 
度 会 进一步 加 剧 。 为 了 实现 土地 利用 的 可 持续 发 
展 ,未 来 土地 管理 决策 必须 对 社会 .经 济 和 气候 因 
素 综合 考量 ,对 祁连山 地 区 不 同 发 展 情景 下 陆地 碳 
储量 进行 预测 ,这 对 祁连山 地 区 碳 储量 增加 和 环境 
的 保护 具有 重要 现实 意义 。 
3.2 碳 储量 经 济 价值 分 析 及 意义 

近年 来 ,全 球 气 候 变 化 愈加 严峻 和 碳 交 易 市 场 
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的 逐步 发 展 ,对 陆地 生态 系统 的 碳 储 量 不 单方 面 注 
重生 态 价 值 ,对 经 济 价值 估算 的 重要 性 日 益 凸 显 。 
中 国 拥有 巨大 的 碳 排放 资源 ,然而 碳 交 易 市 场 起 步 
较 晚 ,发 展 缓慢 ,未 形成 成 熟 的 交易 体系 。 当 前 ,对 
于 碳 价格 还 尚未 形成 统一 的 标准 ,不 同 的 方法 估算 
得 到 的 碳 储 量 经 济 价值 差异 性 较 大 。 例 如 , 胡 峻 周 
等 "9 采用 支付 意愿 法 估算 了 四 川 省 2014 年 森林 植 
被 的 碳 储量 经 济 价值 为 865.75x10 元 , 若 采用 造林 
成 本 法 估算 碳 储量 经 济 价值 为 1517.55x10' 元 ; 张 确 
等 中 采用 造林 成 本 法 估算 了 昆明 山西 森林 公园 碳 
汇 经 济 价值 为 423.6x10' 元 , 若 采 用 瑞典 碳 税率 估算 
碳 汇 的 总 价值 为 1591.10x10 元 。 与 以 往 的 研究 不 
同 ,本 研究 首次 用 复 利 现 值 法 ,估算 了 祁连山 地 区 
的 碳 储量 经 济 价值 。 

在 2000 一 2020 年 ,祁连山 地 区 碳 储量 经 济 价值 
与 碳 储量 均 呈 增加 趋势 ,结论 与 余 娜 等 ”研究 结果 
一 致 ,这 侧面 印证 了 当地 生态 系统 的 生态 功能 和 环 
境 服务 能 力 的 提升 。 碳 储量 经 济 价值 是 由 碳 储量 
和 固 碳 价格 共同 决定 ,同一 时 期 碳 价格 相同 ,所 以 
碳 储量 是 碳 经 济 价值 的 关键 因素 。 但 需要 注意 
的 是 , 碳 储 量 在 2005 年 和 2015 年 出 现 小 幅度 的 下 
降 , 而 碳 储 量 经 济 价值 仍然 呈现 增长 趋势 ,这 主要 
归 因 于 两 个 方面 :第 一 ,虽然 碳 储量 下 降 , 但 由 于 碳 
价格 的 上 涨 ,单位 碳 的 价值 提高 ,使 得 碳 储 量 经 济 
价值 仍然 上 涨 。 第 二 , 碳 价 格 的 上 涨幅 度 远 远 高 于 
碳 储量 的 下 降 率 。 这 表明 ,在 一 定 条 件 下 , 碳 储量 
减少 不 一 定 导致 碳 储 量 经 济 价值 的 下 降 ,可 能 是 因 
为 碳 价格 和 贴现 率 等 因素 的 综合 影响 使 碳 储量 经 
济 价值 上 升 。 

碳 储量 经 济 价值 在 环境 和 经 济 领域 扮演 着 重 
要 的 角色 ,其 为 评估 碳 资 源 的 潜力 ,推动 碳 减 排行 
动 .促进 碳 市 场 的 发 展 以 及 制定 政策 提供 了 基础 和 
参考 中。 未 来 应 当 制 定 并 实施 碳 定价 机 制 . 加 大 对 
低 碳 技术 的 研发 和 推广 .建立 各 利益 相关 方 之 间 的 
伙伴 关系 和 合作 ,完善 全 球 碳 市 场 交 易 。 同 时 政府 
可 以 通过 采取 措施 保护 和 恢复 森林 、 湿 地 和 海洋 等 
碳 密集 型 生态 系统 ,加 强 土 地 管理 和 保护 ,以 确保 
土壤 中 的 有 机 碳 得 到 妥善 保护 和 管理 。 


4 结 论 

选取 2000—2020 年 祁连山 地 区 期 土地 利用 
数据 ,利用 土地 利用 转移 矩阵 和 土地 利用 动态 度 方 
法 分 析 了 土地 利用 变化 特征 ,应 用 InVEST 模 型 估算 
祁连山 地 区 2000 一 2020 年 陆地 生态 系统 碳 储量 , 结 


研 
合 复 利 现 值 公式 估算 碳 储 量 经 济 价值 变化 趋势 。 
得 出 如 下 结论 : 

(1) 2000 一 2020 年 草地 和 未 利用 地 面积 占 比 为 
80% ,是 祁连山 地 区 的 主要 土地 利用 类 型 ,上 且 两 者 转 
化 较为 显著 。 因 城镇 化 的 推进 和 研究 区 特殊 的 地 
理 特征 ,使 祁连山 地 区 在 研究 期 间 土 地 利用 类 型 呈 
耕地 .草地 ,水域 和 建筑 用 地 面积 扩张 ,而 林地 和 未 
利用 地 面积 缩减 的 变化 特征 。 

(2) 2000—2020 年 祁连山 地 区 碳 储 量 增 加 
44.3x10' t, SHR ARH 2% , 呈 "“ 西 北 低 - 东 南 高 ”的 空 
间 分 布 特征 , 延 东南 走向 逐渐 增 大 。 这 归 因 于 人 研究 
区 适宜 的 气候 特征 和 产业 结构 ,使 草地 面积 大 量 增 
加 成 为 碳 储量 增加 的 主要 原因 。 其 中 ,耕地 .草地 
和 水 域 共 增加 73.1x10'*t, 林地 和 未 利用 地 共 减 少 
28.8x10't。 

(3) 2000 一 2020 年 祁连山 地 区 碳 储量 经 济 价值 
增加 1582.9x10' 元 ,增长 率 为 68.1%。 各 土地 利用 类 
型 碳 储量 经 济 价值 占 比 分 别 为 :耕地 3.1% 林地 
19.9% .草地 51.9% , KIN 6.4% 和 未 利用 地 18.6% 。 
其 中 ,草地 的 碳 储量 经 济 价 值 最 大 ,耕地 的 碳 储量 
经 济 价值 最 小 ,草地 碳 储量 经 济 价值 增加 是 祁连山 
地 区 碳 经 济 价值 增加 的 主要 驱动 力 。 
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Impacts of land use change on carbon storage and estimation 


of economic value in Qilian Mountain region 
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Abstract: Human activities are the main drivers of land use change, Which indirectly affects the value of 
ecological services by influencing the structure and function of terrestrial ecosystems. In this article, We analyze 
the characteristics of land use change based on a land use transfer matrix and dynamic attitude method. We 
applied the INVEST model to estimate carbon storage, and we estimated the economic value of carbon storage in 
each period by combining this value with the compound present value formulas from 2000 to 2020 in the Qilian 
Mountain region. The results showed that, between 2000 and 2020, grassland and unused land were the main land 
use types in the Qilian Mountains, accounting for more than 80% of the total area of land, and that the conversion 
these two land use type was the most obvious conversion. The area of forested and unused land decreased 
between 2000 and 2020, and the area of other land types increased in the Qilian Mountains due to the 
intensification of human activities and increased urbanization. From 2000 to 2020, carbon storage was 
significantly correlated with the distribution of land-use types and showed a spatial distribution of “low in the 
north- west and high in the south- east”. In addition, carbon storage increased by 44.26 x 10° t in the Qilian 
Mountains due to the appropriate geographic environment, the strengthening of vegetation protection, and the 
efforts of artificial management. the carbon storage in arable land, grassland, and waters increased by 73.1 10° t, 
while the carbon storage of forest land and unused land decreased by 28.8x10* t. The economic value of carbon 
storage in the Qilian Mountains from 2000 to 2020 increased from 2325.9 10° yuan to 3908.8 10° yuan, with an 
increment of 1582.9x 10° yuan and a growth rate of 68.1%. In particular, the economic value of grassland carbon 
storage increased by 851.8x10° yuan over the 20-year period, a growth rate of 71.5%, which was the main driver 
of the increase in the economic value of carbon storage in the region. Carbon storage is increasing in the Qilian 
Mountains, and the relevant authorities should continue the previous management schemes while remaining alert 
to the impact of the climate change on carbon storage in the future. This study provides a scientific basis for 
setting regional targets for optimizing land use, promoting sustainable development and responding to global 
climate change. 

Keywords: land utilization; InVEST model; carbon storage; economic value of carbon storage; Qilian 
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